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5^ ' Abstract 
Oh! 

<^ ' For some choices of soft SUSY-breaking parameters, the LSP is a stable neutrahno x^i the NLSP 

ly-^ ■ is a chargino Xi almost degenerate in mass with the LSP (Am^ = m~± — m~o ~ m^— few GeV), and 

CN . all other sparticles are relatively heavy. We discuss the potential of a ^/s ~ 600 GeV e^e~ collider for 

. studying such models. 
> ■ 

1. Introduction. As part of the process of planning for future HEP experimental facilities, it is important 

^—1 ' to evaluate as many motivated scenarios for new physics as possible. Certainly, supersymmetry ranks as 

, one of the most successful models of physics beyond the Standard Model, since it can approximately 

, reproduce Standard Model predictions at low energies, while explaining the hierarchy problem. However, 

(~| the mechanism of supersymmetry breaking is not well understood. In general, the different sources of 

^p^' breaking - gravitational interactions, gauge interactions, the conformal anomaly, etc. - lead to different 

Q ' hierarchies of sparticle masses. Many sources may be present at once, so the true model may be quite 
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complicated. Here, we explore a relatively well-motivated set of models in which the gaugino masses are 
non-universal at the GUT scale and, in particular, are such that there is a quite small mass splitting between 
the lightest chargino (xi) and the lightest neutralino (xi); which is the LSP. Small Am^ = ?7i~± — m~o 
occurs when M2 < Mi <^ which arises naturally when the gaugino masses are dominated by or entirely 
generated by loop corrections. Models of this type include the 0-H superstring model |2|, |3|, ^ and the 
(AMSB) models in which SUSY-breaking arises entirely from the conformal anomaly |5|, ^, 0]. The same 
hierarchy also occurs when SUSY is broken by an F-term that is not an SU(5) singlet but rather is a 
member of the 200 representation contained in (24x24)symmetric = 1 e 24 e 75 e 200 [|. Techniques for 
detecting and studying supersymmetry are very dependent upon Am^ and on the relative magnitude of 
m~± ~ m~o as compared to other supersymmetric particles. Previous discussions of small Am^. scenarios 
at lepton colliders include early studies in the context of the 0-H model ^, Q and various studies 
specific to the AMSB boundary conditions ||9|, |ll|, 12 1. This study will focus on signals and parameter 



determination for small Am^ at a ~ 450 — 600 GeV, L ~ 1 ab^^ linear e+e^ collider (LC) when the 
only SUSY particles kinematically accessible at the LC are the Xi ^iicl Xi^. 

To illustrate the possibilities for Am^ and other gaugino masses, it is useful to give some tree- level 
results for the gaugino mass parameters at the scale mz- The O-II model with 5gs = ~4 yields M3 : 
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M2 : Ml = 6 : 10 : 10.6, the O-II model with 6gs = (equivalent to the simplest version of the conformal 
anomaly, AMSB, approach) yields M3 : M2 : Mi = 3 : 0.3 : 1, while the 200 model yields 6 : 4 : 10. As a 
result: in the 200 model and the 0~II 5gs = (or pure AMSB) model, M2 <C Mi and Am^ = ""T-^i ""^^o 
can be smaller than ~ m.,^ at tree-level; in the O-II 5gs = ~4 case, M2 is only slightly less than Mi implying 
that Amy- < a few GeV is very typical — Am^ < 1 GeV if > 1 TeV as when RGE electroweak symmetry 
breaking is imposed |p. Loop corrections can be significant. In particular, even if Am^ < rn-jr at tree- 
level, radiative corrections ||9|, 14 1 usually increase Am^^ to above m^r (see e.g. the Aniy- graphs in Q); 
Aniy. < niT^ is only possible for very special parameter choices. Most typically, Am^ is predicted to lie in 
the range from slightly above m-,^ to several GeV. 

All other SUSY particles could have masses substantially above m~± . In particular, in the O-II 5gs = 

Xi 

and AMSB models, m~o is significantly larger than m~±. Only in the O-II 6gs = ~4 model is m~o close 

X2 Xi X2 

to In most cases, large |//| is required by radiative symmetry breaking so that the X2, X3 Xa 

are very heavy. The masses for the squarks and sleptons are uncertain, but could be quite large, which 
is the case on which we focus. In particular, we neglect sneutrino exchange contributions to the xtxi 
and "yxtxi cross sections. For example, both cross sections start to become increasingly suppressed for 
y/s ~ 600 GeV as the electron sneutrino mass is decreased below 1 TeV; see, for example, ||l5[| . 

The neutralino and chargino couplings to W and Z bosons in the wino-LSP scenario were reviewed in 
|11|. One finds that the ZxiXij ■^XiX2) ^X2X2^ and VF^xf X2 couplings (and corresponding cross sections) 
are all small, while e~^e~ ^ Z, 7 ^ xtxi^ixtxi and e~^e~ W^xfxi can have large rates. 



Am^{ MeV) 
cr( cm) 


125 
1155 


130 
918.4 


135 
754.1 


138 
671.5 


140 
317.2 


142.5 
23.97 


150 
10.89 


Am^{ MeV) 
cr( cm) 


160 
6.865 


180 
3.719 


185 
3.291 


200 
2.381 


250 
1.042 


300 
0.5561 


500 
0.1023 


Am^( MeV) 
cr(cm) 


600 
0.055 


700 
0.033 


800 
0.019 


900 
0.011 


1000 
0.0072 


1500 
0.0013 


2000 
0.00036 



Table 1: Summary of cr values as a function of Am^ 



The most critical ingredients in the phenomenology of such models are the lifetime and the decay modes 
of the Xi^, which in turn depend almost entirely on Am^ when the latter is small. The cr and branching 
ratios of the Xi^ as a function of Am^ have been computed in Q. For Am^ < mj^, only Xi ~^ e^z^eXi 
is important and cr > 10 m. Once Am^. > myr, the x^ ~^ '''"^Xi mode turns on and is dominant for 
Arrij^ < 800 MeV, at which point the multi-pion modes start to become important: correspondingly, one 
finds cr < 10 — 20 cm for Am^. just above 771,^ decreasing to cr < 100 /um by Arriy. ~ 1 GeV. For later 
reference, we give some specific values of cr as a function of Am^^ in Table |l|. 

Finally, we wish to comment on what values for m~± are most natural in the context we are considering. 

Xi 

The constraint from fine-tuning, though somewhat ill-defined, does provide some guidance. The degree of 
fine-tuning is largely controlled by the magnitude of the gluino mass. In the models discussed above there 
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are considerable differences in tlie ratio of Mj^jM^- This ratio is largest (~ 9) in the {bcs = 0)/AMSB 
scenario and fine tuning increases rapidly with M2 ~ ni-i to levels that are highly problematical once 
m~± > 200 GeV. In the other scenarios M3 is close to M2 and fine-tuning is a less severe concern. The 

Xi 

main focus of this paper will be on scenarios in which the chargino and neutralino are highly degenerate, 
which is most naturally the case in the {bcs = 0)/AMSB scenario. Thus, our focus will be on procedures 
relevant for m~± < 200 GeV and Amy < 1 GeV. 

Xi ^ 

2. General Discussion of Detector and Signals. For our discussions, we consider a detector with the 
components listed in Table (See [0] for further details.) The SVX, CT and PS all give (independent) 
measurements of the dE/dx from ionization of a track passing through them. This makes it possible to 
distinguish a heavily-ionizing chargino (which would be > twice minimal ionizing [2MIP] for (3^ < 0.85) 
from an isolated minimally ionizing particle [IMIP]. The net discrimination factor would probably be of 
order few x 10~^. In our simulations, we employed an efficiency of 90% for tracks with 13^ < 0.85 p7| ]. 

The possible signals based upon detecting a non-promptly decaying xf were detailed for a hadron 
collider in |Tl|. Appropriately modified versions for XiXi 'yxtxi production at the LC are listed in 
Table |. 



Component 


Description 


SVX 


Silicon vertex detector from close to beam pipe to ~11 cm. 


CT 


Central tracker starting just past SVX. 


PS 


Pre-shower just outside the tracker. 


EC 


Electromagnetic calorimeter. 


HC 


Hadronic calorimeter. 


TOF 


Time-of-flight measurement after HC and just before MC. 


MC 


Muon chamber with first layer after the HC and just beyond the TOF. 



Table 2: Summary of detector components referred to in the text. 



Assuming that only the Xi Xi light, SUSY particle production will be primarily in the final 
states xtxi J IXiXi 7 and (when not phase space limited) Vl^^XiXi^- (The 7X1X1 cross section is suppressed 
due to the small ZxiXi coupling.) The possibilities for SUSY detection depend largely upon which aspects 
(if any) of a Xi^ in the final state are visible ||2|, |3|, ^, which in turn depends almost entirely on Am-j^. We 
distinguish several interesting ranges for Am^: 

• Aniy. < myr: the xf yields a 'stable particle' LHIT and/or DIT track and is easily detected: xtxi 
production will be easily seen. 

• ttItt < Amj^ < 1 GeV: the Xi^ — > Xi^^ decay yields a soft vr track, possibly in association with a 
STUB [Aniy, < 180 MeV) or HIP (Am^ < 1 GeV) signature. Direct Xi^xr production then yields a 
^ + TTvr final state, where the is associated with the Xi's- Since the vr's are very soft, backgrounds 
to this final state from 77-induced interactions are very large, and it is unlikely that the tt^tt' SUSY 
signal can be isolated. One must tag Xi^xr production using e"'"e~ 7Xi^xr' leading to a "f + I^ + tttt 
final state, or employ the more kinematically limited W^Xixt ~^ (^^) ?) + + '^^ final states. 
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Signal 


Definition 


LHIT 


At least one long, heavily-ionizing (> 2MIP's as measured by SVX+CT+PS), large- 
Pt track that reaches the MC. The energy deposit in the HC in the track direction 
must be consistent with expected ionization energy deposit for the /? measured (using 
TOF and/or SVX+CT+PS), i.e. no hadronic energy deposit. 


TOF 


At least one large-pr track seen in the SVX and CT along with a signal in the TOF 
delayed by 500 ps or more (vs. a particle with /3 = 1). HC energy deposit (in the 
direction of the track) is required to be consistent with the ionization expected for 
the measured /?, i.e. no large hadronic deposit. 


DIT 


At least one isolated, large-pr track in the SVX and CT that fails to reach the 
MC and deposits energy in the HC no larger than that consistent with ionization 
energy deposits for the measured (using SVX+CT+PS) (3. Heavy ionization in 
the SVX+CT+PS, corresponding to /? < 0.8 or (3 < 0.6 (DIT8 or DIT6), may be 
required. Large ^. 


KINK 


At least one track that terminates in the CT, turning into a soft, but visible, 
charged-pion daughter-track at a substantial angle to parent. Large 1^. 


STUB 


At least one isolated, large-pT (as measured using SVX) track that registers in all 
SVX layers, but does not pass all the way through the CT. Energy deposits in the 
EC and HC in the direction of the track should be minimal. Large I^. 


SNT 


One or more STUB tracks and large 1^ with no additional trigger. Heavy ionization 
of the STUB in the SVX corresponding to (3 < 0.8 (SNT8) is required. 


S7 


One or more STUB tracks and large ]^ with a pj. > 10 GeV trigger. Heavy ionization 
of the STUB in the SVX corresponding to (3 < 0.8 (S78) is required. 


HIP 


At least one high-impact-parameter {b > 5ab) soft pion track in the SVX, with 
> 10 GeV triggering and large perhaps in association with a visible KINK in 
the SVX. 




Isolated, large-ii^T photon and large I^. Relevant if the soft vr's cannot be detected. 


mSUGRA-like 


jets or leptons + ]f 



Table 3: Summary of signals for XiXi and 7X1 Xi production. MIP refers to a minimally-ionizing-particle 
such as a highly relativistic muon. For detector component notation, see Table || 
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• 1 GeV < Am^ < 2 GeV: the Xi decays with roughly equal probability to Xi'^^^ Xi'''"^'''"'' and Xi^^'^ 
{£ = e,fi). The pion(s) or charged lepton will be rather soft and a 7 tag will probably still be 
necessary to eliminate backgrounds. 

• Am^ > 2 — 3 GeV: the Xi decays either to x5+™ulti"Pioii modes at the low end which start to 
resolve into jets at the higher end or to Xi^^^- The visible decay products are sufficiently energetic 
that 77 induced backgrounds can be rejected (using a combination of event topology and ]^) to the 
extent necessary for mSUGRA-like mode detection of direct e~^e~ XiXi production. 

Sensitivity of LEP2 detectors to these various signatures was sufficient to exclude [|^, [l^: (a) m~± < ■v/s/2 
in the 'stable' and 'standard' regions of Am^^; and (b) w-^i < 80 GeV (assuming the z^e is heavy) in the 
"T-TT < Am^ < 2 GeV region. In particular, the backgrounds to the 'stable' and 7-tag+^+soft-7r signals are 
very small. In our analysis, we assume that they can continue to be neglected relative to our signal rates at 
the LC after simple cuts. For the 7-tag+soft-7r signal, this may, at ffi'st sight, seem problematical given that 
the signal cross section declines with increasing s, whereas the two-photon (7*7*) collision backgrounds 
tend to increase logarithmically with s. Below, we review the simple cuts that are likely to be very effective 
in eliminating the two-photon backgrounds. 

It is first important to understand the kinematics of the signal. We will trigger on a photon with 
Pj^ > 10 GeV and 10° < 6.y < 170°, 10° being the angle at which electromagnetic coverage by the detector 
is expected to begin. We will require that = (pg+ + Pg- —p^Y, the invariant mass-squared of the 
XiXi pair, be such that M^^ > 2m~± when searching for Xi^'s of a certain mass. Of course, the search 
will start at fn^± > 80 GeV, the LEP2 limit. Most of will be invisible, being carried by the XiXi 
pair. The 7r+ and vr" in the final state will be very soft, with energies basically set by the size of Am^^, 
and largely central, with some bias for the vr"^ (tt") to move in the direction of the e~ (e^) beam. Thus, 
they will tend to be somewhat back-to-back. However, their angle of acoplanarity will be very uniformly 
distributed. Finally, the pTp of the trigger photon will be clustered at the lowest allowed values and will be 
primarily balanced by the missing transverse momentum, carried by the Xi's- 

There are two different types of two-photon background to the 7 + vr+vr^ + Ip final state of interest. 
First, there are e^e~ — > 7 + e^e~ + 7*7* ^ 7 + e'^e~ + X, with X = vr^vr^, reactions in which the 
final e"*" and e~ are both lost down the beam pipes and, thereby, provide the large missing energy and 
missing mass associated with the XiXi pair in the signal reaction. If we trigger on (i.e. tag) a photon with 
substantial pj., either the e~ or the will be given a sufficient transverse kick that it will be detected. 
To be precise, we define 9^ to be the angle above which an electron or positron can be detected. The 
largest transverse momentum, p™^^, that can be carried by the final e^e~ pair without one being detected 
arises when the final and e~ are coplanar, have relative = and have polar angles with respect 
to their respective beams, 9^ and 0~ , both equal to 9^: p^^^ = [E^ + E^) svnOd. Using the constraint 
E"^ = ^/s - E~ - E"^ - E^ , we have p^'^^ = (^s - E^ - E'^) sinOd- Meanwhile, the smallest pr of the 
the 7 + X system in the final state is p™™ = Pt ~ Pt ' where denotes the magnitude of the X system 
transverse momentum, p™™ < p!^^^ is required to avoid detection of both the final and e~ . This 
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condition can be rewritten in the form 

< • 

For the scenario under consideration, and are of order ISm^ times a modest boost factor, and do 
not exceed 1 GeV for the scenarios of interest. Thus, if we impose a pji cut that exceeds -v/sj^^^ then 
this background process cannot contribute. Current detector plans are consistent with instrumenting the 
beam hole down to ~ 1°, for which a cut of pj, > 0.0172-y/s will suffice. For ^fs = 600 GeV, we would 
require > 10.3 GeV, which motivates our nominal minimal cut of > 10 GeV. With this cut, the only 
way in which an e~^e~ ^ 7 + e'^e" + X process can contribute to the 7 + tt^tt" + ]f final state of interest 
is if X itself contains missing energy of magnitude approaching the magnitude of the minimum being 
required. One candidate is X = t~^t~ 7r~^TT~iyV, where the t~^t~ invariant mass is large enough that 
the mass and momentum of the lyV system is at least 2 to 3 GeV. However, the vr^ and tt~ will tend to be 
much more energetic than expected in the signal reaction (unless Am^ > 500 MeV) and will tend to be 
produced more forward and backward (with correlations with the beam directions opposite that expected 
for the signal) because of the i^-channel r exchange in the Feynman diagram that gives rise to 7*7* — > t~^t~ . 
Further, the background from this process can be independently measured using the t ^ pv decay modes, 
and then subtracted. Another candidate is X = tt^tt^ Z{-^ uV) for which e^e~ ^ 7 + e^e~ + X has a 
small cross section. In events of this type, virtualities are large and the vr's produced will tend to be very 
energetic and largely forward and backward, whereas the vr's from the signal reaction are quite central and 
quite soft in the lab frame. Cuts that require soft central pions will thus be very effective in eliminating 
such backgrounds while retaining almost all the signal events. 

A second type of two-photon background arises from overlapping e^e~ ^7^ (e.g. ^Z*{-^ I'V), with 
mz* > ^xx) ^'^^ e+e^ e+e^7*7* e~^e~TT~^TT~ collisions in which the final e^e^ of the 2nd reaction 
are allowed to disappear down the beam lines regardless of pj.. At TESLA, with large bunch separation, 
the two e^'s must come from the same bunch and the two e+'s from the same e"*" bunch. Note that 
a cut requiring large I/!, in particular large ^xx' only be satisfied if the ^ comes from the 

same reaction as the 7. The 'hard' process will have a small cross section. Further, the vr's from the 
overlapping two-photon event tend to be very energetic and the cut requiring central/soft vr's will be very 
effective (and essentially 100% efficient for the signal events). Further, the tt"*" and tt" produced in the 2nd 
collision will tend to balance each other in transverse momentum, whereas this will never be the case for 
the vr's from the signal events. In addition, for most such events there would be substantial separation in 
z (displacement along the beam line) between the two overlapping events. The LC vertex detectors will 
have exceptional ability to determine the z vertex location of the event producing the tt's, assuming cr 
for the xt decay is not large. (If it is large, observation of large transverse impact parameters for the tt's 
would, alone, be sufficient to eliminate any background.) The z location of the collision producing the 7 
will hinge on electromagnetic calorimetry and will require a pre-shower plane for good accuracy. We are 
uncertain of what resolution to expect here, but it could be sufficiently good that a non-zero z separation 
would be seen for most overlapping collisions, which could be used to eliminate most of this background. 
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Another cut that is useful[| for ehminating two photon backgrounds is to require that the missing energy 
vector ^ = pg+ + pg- — P'y — ~ Pn- have an angle of > 100 mrad with respect to both beam axes. In 
fact, our cuts of pj^ > 10 GeV and 10° < 9j < 170° guarantee that this is the case for all signal events; 
indeed, all but a tiny fraction of the events have an angle between ^ and both beam axes of > 200 mrad. 

In the absence of a reliable Monte Carlo for generating these types of backgrounds, it is impossible to 
be absolutely certain that a collection of small contributions will not add up to a non-negligible result. 
The true background level will probably only be known once data taking has begun. If a background is 
found for our minimal tag requirements of pj. > 10 GeV, 10° < 9^ < 170° balanced by large missing mass 
and soft central pions, the pj- cut, can, for m~± not too near \^/2, be increased to the point where the 
background disappears. For L = 1 ah^^ of integrated luminosity, this will typically have little impact 
on our ability to discover 'yxtxi production and determine m~± and Am^^. In fact, one very good set 
of procedures for the latter determinations, discussed in a later section, focuses on events for which the 
photon energy is near maximal and the invariant mass of the xtxi system is very close to 2m~±. In this 
kinematic region, backgrounds will surely be negligible unless m~± is close to \fs/2 (implying that there is 

Xi 

no kinematic room for increasing the photon tag requirements). To be more quantitative, we show in Fig. ^ 
the unpolarized beam and right-handed polarized e~ beam cross sections as a function of the minimum pj, 
accepted. Results are shown for the representative cases of m~± = 175, 250, 290 GeV and n = 300, 600 GeV. 

Xi 

In the unpolarized case, L = 1 ah~^ gives acceptable event rates for discovery (we assume 10 events will be 
adequate if there is no background) even for very large pj, cuts so long as m~± < 275 GeV. As discussed 
in a later section, these event rates for large pj, are also quite adequate and, indeed, very useful for precise 
determinations of and Am^ for such xf masses. (The ability to determine ?7i^± from the event 
rate as a function of the pj, cut is immediately apparent from the figure.) At 'm~± = 290 GeV, a cut 
of pj^ > 17 GeV would still leave a (discoverable rate) of 30 events (for L = 1 ab^^) if the background 
is negligible. Determinations of ^n-ji and Arrij^ would become less accurate at such a high mass. In the 
polarized case, event rates are very sensitive to the n parameter. This is good in that it provides important 
sensitivity to the (and M2) SUSY parameters, as described in a later section. However, Fig. |l| makes clear 
that the polarized cross sections are not necessarily large enough to be accurately measured and might not 
be observable at all. Use of the polarized cross sections will, in many cases, require that the background 
really is very small starting at pj, ~ 10 — 20 GeV. Finally, we note that the results of Fig. || assume that 
the amplitude contribution from the sneutrino exchange diagram is small. If the sneutrino mass is at or 
below ^/s, there will be substantial suppression of the cross section due to the z/g exchange diagram. 

We reemphasize the fact, apparent explicitly in Fig. |2| (discussed in the following section), that there 
is a high probability that one or both of the pions in the "yxtxi final state will have an observable impact 
parameter whenever Amy is below ~ 1 GeV and m~± is in the m~± < 200 GeV region, as preferred for the 

^ Xi Xi 

i^GS = 0-II)/AMSB boundary conditions. In this case, backgrounds will unquestionably be negligible 
and all the parameter determination studies described later will be possible with the described precisions. 
^We thank H.-U. Martyn for pointing this out. 
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Vs=600 GeV, a Dependence on p| Cut 




Figure 1: The ^/s = 600 GeV cross sections for JXiXi production as a function of the minimum photon 
Pj^. The results for an unpolarized and a right-handed polarized e~ beam are shown in the upper and 
lower plots, respectively, for the cases of m~± = 175, 250, 290 GeV and fi = 300, 600 GcV, with tan/3 = 5. 
Cross sections arc computed requiring 10° < 9^ < 170°. These results assume that the sneutrino exchange 
diagram can be neglected. 

Finally, we return to the possibility of using e^e^ ^^XiXi which the W can decay either 
leptonically or hadronically and xt ~^ ^^Xi- Fo'^ moderate m~±, the cross section for the production 
process is quite substantial. For example, for m~± = 175 GeV, tan /? = 5, = 600 GeV and ^/s = 600 GeV 

Xi 

(450 GeV), the unpolarized beam cross section is 12 fb (0.46 fb) and the pure cross section is 0.28 fb 
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(0.007 fb), after summing over both W~XiXi and W~^XiXi production. The ^/s = 600 GeV cross sections 
(for instance) are to be compared to the pj^ > 10 GeV cross sections of Fig. [l|: 38.5 fb for unpolarized 
beams and 0.27 fb for pure e^. The W^xtxl cross section is thus substantial so long as 2m~± + niw is 
not too near ^/s. This production mode could certainly be considered as an alternative for discovery and 
parameter studies if the "yxtxi mode should (contrary to our expectation) have significant background. 
We have not studied possible backgrounds to the W^xfxi channels, although we expect backgrounds to 
be small. We will only remark on its possible utility in a few places. 

Vs=600 GeV 



standard SUSY Search 50 fb~^ 






100 120 140 160 180 200 220 240 260 280 300 

Chargino Mass (GeV) 



Figure 2: Regions in Amy, m~± parameter space for which different discovery modes are viable. The SNT, 
7+HIP and 7 + vrvr modes are presumed to be background free; we assume that 10 events will be adequate 
for discovery. For the SNT, 7+HIP and 7 + signals, discovery reach is given for L = 50 fb~^ and also 
for L = 1 ab~^. Also shown by the vertical band and vertical line (the band for L = 50 fb~^ and the 
line for L = 1 ab~^) is the reach of the 7 + detection mode, which is relevant only if the soft vr's are 
not detectable, as might be the case for some range of Am^ ^ 200 MeV. (Limits in this mode assume 
S/^/~B > 5 and S/B > 0.02 is required; see text.) Note the large increase in discovery reach for the 7 + 
mode with increased L. 

3. Discovery Reach. Following the above discussion, we can estimate the regions in (Am£,m~±) 
parameter space for which the various SUSY discovery modes will be viable at a y/s = 600 GeV LC. They 
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are shown in Fig. ||. For Aniy. > 2 GeV, we assume (fohowing the LEP experience) that observation 
of XiXi production will be possible up to very nearly n^^± ^ ^/sf^- For Am^ < mj^, the xf will be 
sufficiently stable in the detector that an LHIT signal for XiXi production will be easily detected up to 
threshold. For Arriy. between and roughly 200 MeV, the SNT signal will be viable. In fact, since 
triggering will not be necessary at the NLC, even (STUB) charged tracks as short as 4 or 5 cm can be 
directly imaged, f\ implying that our /? < 0.8 heavy ionization requirement for the SNT signal could be 
relaxed. For 200 MeV < Am^ < 2 GeV, it is very probable that one will need to employ the 7 tag 
signatures, i.e. 7+HIP(s) and 7 + 7r(s) (always with large 1^). We believe that the 7+HIP and SNT 
signatures will be background-free. (They disappear at the largest values of m~± due to inadequate boost 
for the produced Xi, which leads to too small an impact parameter for the HIP or too short a path 
length for the STUB track.) Development of Monte Carlo programs that can accurately compute the 
backgrounds to the above signals should be a high priority, especially for the 7 + ^ + 7r(s) (7 + 7r(s) for 
short) signal. However, the measured background after cuts at LEP2 for the 7 + 7r(s) signal was negligible 
or very small even without requiring a high impact parameter for at least one of the tt's. This, coupled 
with the arguments given in the previous section allow for some optimism that the background for the 
7 + 7r(s) signal will continue to be small at high ^/s, even for the relatively mild pj^ > 10 GeV cut. The 
exact boundaries for the LHIT, SNT, 7-|-HIP(s) and 7 + 7r(s) signatures shown assume that 10 events are 
adequate to detect the above signals (assuming no background). The rapid s-wave threshold turn-on of 
1X1X1 production allows one to probe almost to the fn~± ~ ^fsjl — Pj, kinematic limit so long as the 
detector can resolve the soft vr's from the Xi decays. As shown in Fig. ^, even assuming zero background, 
L = 1 ab^^ does not greatly increase the discovery reach of the STUB and HIP channels as compared to 
L = 50 fb-^ The increase in reach for L = 1 ab ^ for the (stable) LHIT signal is very small and is not 
shown. Since it is possible that backgrounds that are negligible at L = 50 fb~^ result in a non-negligible 
number of events at L = 1 ab~^, one should not absolutely rely on the indicated increase in coverage for 
the SNT and 7-hHIP signals for the latter L. 

Given that the boundary conditions in models with loop-dominated gaugino masses are such that 
Arriy. £ [200 MeV, 2 GeV] is almost a certainty, the reach of the 7-tag-l- ^-l-soft-Tr's (with or without 
HIP(s)) signal is of great importance. If the soft pions cannot be detected/tracked before they curl up (e.g. 
because the magnetic field is too strong), then sensitivity to jXiXi production will be much more limited. 
(This is illustrated by the somewhat higher m~± limits of L3 vs. DELPHI for Am^ e [200 MeV, 2 GeV] 
due to the fact that L3's magnetic field is weaker than DELPHI'S.) If the vr^ tracks are not detected, one 
important signal is e"*'e~ IXiXi ^ 7 + invisible, for which the e+e~ ^vV background is substantial. 
For any given presumed value of m~±, we place a cut on the invisible mass, 4^, recoiling opposite the 
trigger photon of at least I^I > 2m~±, that removes much of the z/F background, in particular that from 

Xi 

on-shell ^ + Z production. The more limited discovery reach is illustrated by the 7 -|- iyf boundaries (one 
for L = 50 fb~^ and one for L = 1 ab^^) drawn at fixed for Am-^ G [200 MeV, 2 GeV] in Fig. For 

■^We thank M. Peskin for bringing this to our attention. 
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Figure 3: dN/dl^I for signal and background in the 7 + i^f final state using 5 GeV bins and assuming 
L = 50 fb^^. Also shown are S and 



lower m~± values, cuts can be placed on tfi so as to obtain S/B > 0.02 and S/^fB > 5 for L = 50 fb ^. 
The difficulty of observing the signal for large m~± can be seen from the m~± = 200 GeV, L = 50 fb^^ 

Xi Xi 

results shown in Fig. ^. The upper part of the figure gives S and 5 + S in 5 GeV bins in the > 525 GeV 
region. The lower part of the figure shows S and ^/B for 5 GeV bins. We see that S/ \fB per bin is only 
of order 1 for 4// > 525 GeV. Although at the peak in S", located at 4^ ~ 580 GeV, the 5 GeV bin has 
S/B = 0.028, for L = 50 fb^^ this one bin only yields S/VB = 1.87. Thus, for this L it is necessary to 
include all the bins with > 525 GeV in order the get a net S/^/B > 5 signal; in this case the net S/B is 
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~ 0.02. However, with higher L one can do much better by focusing on the highest bins for which S/B 
is largest. In Table we give the 7 + 4^ cross section cr as a function of m~± for 580 GeV < l/I < 590 GeV, 



m~± (GeV) 

Xi 


175 


200 


250 


260 


275 


285 


290 


a(fb) 
S/B 
L (fb-i) for 5fT 


3.9 
0.043 
150 


3.6 
0.040 
175 


2.3 
0.025 
435 


1.8 
0.020 
690 


1.4 
0.016 
1140 


0.90 
0.010 
~ 3000 


0.52 
0.006 
~ 8500 



Table 4: 7 + 4^ results for 580 < 1\/I < 590 GeV. 

the ratio S/B^ and the luminosity required for S/\/^ = 5. From this table, we see that if S/B > 0.02 is 
required for a viable signal, then m~± up to ~ 260 GeV can be probed for L > 700 fb~^. If systematics 

Xi 

can be controlled to the S/B 0.01 level, then the 7 + iyf signal could be observed up to m~± ~ 285 GeV 

Xi 

for very high L ~ 3000 fb^^. 

4. Determining m~± and Am^. If the chargino or its decay products are directly visible, for expected 
e+e^ collider luminosities it will be possible to detect xtxi or 'yxtxi production for m~± up to values 

Xi 

very close to the \fs/2 threshold. The soft 7r(s) from the decay or the non-prompt xf decay LHIT, 
SNT and/or 7-I-HIP signals will indicate clearly that Am^ is small. The next important task will be to 
measure m~± and Amy as precisely as possible. We focus entirely on scenarios with 771,7^ < Amy < 0.8 GeV 



for which we must employ the "yXiXi channel and for which the Xi ~^ '''"^Xi branching ratio is close to 
1. For larger Am,^, one would have to take into account the Xi^ — > vr^vr'^x^ xt ~^ ^^^Xi decays. Our 
techniques could still be applied, but the effective event rates for the ir^Xi channels would be reduced and 
some accuracy lost unless the other modes could be included in the analysis following analogous procedures 
(which are relatively straightforward extensions of those presented). 

The most direct technique for determining it^:^± and Am^ independently of theoretical assumptions 
is via the vr^ kinematic distributions. We shall later discuss the determination of Am-j^ via observing a 
distribution of vr^ impact parameters. By way of preview, we note that the impact parameter distribution 
is determined primarily by cr and Amy, with weak dependence on m~±. Once m~± is determined via the 

Xi Xi 

kinematic procedures we shall describe below, in favorable cases (large cr) the impact parameter distribu- 
tion will allow a simultaneous determination of cr and Am,^ with errors on Am,^ that are competitive with 
the errors we shall find using kinematic distributions. Of course, the kinematic distribution procedures we 
now describe are essential when the Xi decay is simply too prompt to directly observe its path length. 
However, even if the distribution of path lengths is not very apparent (because most events are clustered at 
low impact parameter, e.g. as typical for Am^ ~ 800 MeV for which cr ~ 190 fim) or useful, most events 
will have an observable non-zero impact parameter for at least one vr and, consequently, backgrounds to 
the 7 -|- TT^-K~ + }/! signal will certainly be negligible after our simple cuts on pj,, 6^ and M^^ described 
earlier. 

The observables in the "fxtxi ~^ 7^^"^ XiXi final state of interest are the invariant mass of the 
chargino pair, M^^ = J {Pe+ + Pe- "^7)^1 and the energies, E*, of the vr's in the M-^^ rest frame. For 
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LC with Vs = 600 GeV 




350 400 450 500 550 600 

Figure 4: We plot event number contours in [M^-^, E*] parameter space for several different values of Am^^ 
taking Tn-i = 150 GeV. Note the event concentration at high M^^. The arrows indicate the values of 
Am^^, which are approximately equal to £'* at the M^^ = 2m~± endpoint. The vertical line at 360 GeV 
denotes the M^-^ value below which the E* distribution is approximately Gaussian. 



It 

Xi rest frame with respect to the boost direction. The maximum and minimum E* values are given by 



M^y^ > 2m~± , the Xi are boosted in the M^^ rest frame and E* depends upon the orientation of p-,^ in the 



£;j^-- = 7~±(i?o + /3~±Po), El-^'^ = -f~±{E,-l3~±p,), (2) 
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with 



4m~+ -1 m~+ + ml — m\ 



characterizing the charginos in the M^^ rest frame and the vr^ in the rest frames. These equations can 
be inverted to give 



2 e"^ + ml- E"^ ±J{E^ + ml- E'^Y - AE^ml 
^xf = W ' 



Xi 



2m~± ^ 
m~± h m| ~ , (5) 



27~± ^ Xi \ xi 7~± •■ mr± 



where 2^ = E^ max_^^^ min g^j^^ = E* -^^^-E^ (The ~ expressions in Eqs. (|) and (D apply if we 
neglect terms of order m^/m~± and Am|/m~±.) For any given M-^^, there is a two- fold ambiguity in 7~±. 
The — (+) sign applies for lower (higher) M^^ values. Only one set of choices as a function of will 
lead to M^^-independent values for and Am^. Even at a fixed M^^, fitting the actual distribution in 
E* can provide a unique determination of "t-^± and Am^, given enough statistics. 

Since the vr's are soft, the resolution for measuring p-,^ in the lab frame will be dominated by the 'con- 
stant' term, which has a typical value of ~ 0.5%. Further, bM^^jM^^ ~ bE^^jM'^ ~ O.iyE^-y/i/M? 



(all in GeV) will also be very small. We have generated events including the measurement smearing in E-^. 
For each event with a given M^-^, we boost to the M^^ (i.e. chargino-pair) rest frame and compute the 
energies, -E*, of the observed vr's. In Fig. ^, we show the region of the [£^*,Mj^^] plane occupied by the 
events for several different choices of Am^^, taking W'~± = 175 GeV. For L = 50 fb^^ and = 600 GeV, 
the occupied region will contain ~ 5000 entries (2 entries per event). A large fraction of the entries reside 
near the large-M^^ boundary. 

The location of the threshold in provides a particularly useful way to measure 2m~± . The statistical 
error for such a determination depends upon the number of events near the threshold. We assume that 
the 7 + TTTT signal is background-free and compute the approximate \a error bm~± from the criterion: 

Xi 

S{m~± - bm~±) - S{m~±) 



^S{m~±-bm~±] 



where S is the number of events in the near-threshold region 2m~± < Myy < 2m~± + 10 GeV. For 
m~± G [150,225] GeV (g [250,275] GeV) we find bm~± ~ 1 GeV (~ 0.5 GeV) for L = 50 fb^^ and 

Xi Xi 

6m~± ~ 0.2 GeV (~ 0.1 GeV) for L = 1 ab^^. The accuracy improves at higher masses because of the 

Xi 

increasing sharpness of da/dMy.y. as a function of M^^. 

The most direct way to determine Am^ is to concentrate on the ~ ^"^x"*" '^^Sion for which 

E* ~ i?o — Am^. Because of resolution smearing, there are actually events with M^^ < 2m~± and a 
finite spread to E* in this region. We have examined the E* distribution for two cuts: (I) < 2m~± 
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Figure 5: We plot the distribution of for the case of m--± 



175 GeV, Am^ 



200 MeV and M^^ < 



360 GeV. Also shown is the Gaussian fit and its r.m.s. used to determine the error for this case in Table |. 
We took L = 1 ab~^, yielding 392 events in the plot. 



and (II) M^^ < '^rn~± + 10 GeV. For cuts (I), there are fewer events but, for Am^ < 400 MeV, the E* 
distribution is more closely centered on Am-^. A typical case is illustrated in Fig. |5|, based on cuts (II). 
The error on E*, the average vr energy, is given by 6E*/^/N, where N is the number of events and 5E* is 
estimated as the width of a Gaussian fit to the E* spectrum. Table || gives E* and SE* for both M^^ cuts 
(I) and (II) for a number of "i^± and Am,^ choices assuming L = 1 ab~^. L = 50 fb~^ is not adequate to 
give particularly small errors, but would allow a first determination of Am^ to within 10 to 20 MeV or so. 
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cut I 




cut II 


[m~± (GeV), Am^ (MeV) J 


E* (MeV) 


f T — ! J. / TV T T 7'\ 

oE* (MeV) 


7 — 1j- /TV /r T 

E* (MeV) 


C 7 — 1j- /TV /r T 

oE* (MeV) 


[175, 200] 


201.6 


1.9 


203.8 


0.9 


[175, 300] 


303.4 


3.4 


306.7 


1.8 


[175, 500] 


493.2 


9.0 


499.0 


3.6 


[225, 200] 


200.6 


1.3 


203.2 


0.7 


[225,300] 


300.4 


2.6 


303.6 


1.2 


[225, 500] 


493.9 


6.8 


498.3 


2.7 


[250, 200] 


200.4 


1.0 


202.8 


0.5 


[250, 300] 


299.9 


2.1 


303.5 


1.0 


[250, 500] 


494.0 


5.0 


498.2 


2.1 



Table 5: Summaxy of values and errors, SE*, as a function of n^^±, Am^ and M^^ cut, for m~± = 175, 
225 and 250 GeV and Am^ = 200, 300, and 500 MeV, taking L = 1 ab~^ 



There is considerable variation of 5E* with the case and the cut. In addition, E* is usually not precisely 
equal to Am^, with shifts ranging up to 7 MeV, depending on the exact case and cuts employed. However, 
the expected shape for the E* distribution for any given choices of nT-^±, ^'^x ^^'^ ^xx precisely 
known (to the extent that the resolutions for E.,^ and E^ are known) and the expected shift E* — Amy. can 
be computed. Our estimate is that one could in the end achieve an uncertainty for Am^ of order 2-5 MeV 
or better in most cases. Once actual data were available, one would take the experimental distribution 
resulting from the underlying J^i-i and Am^ and compare it to a selection of theoretical predictions for 
the measured m~± (using the technique described earlier) and a range of Am^ choices and minimize the 

To repeat a point from the introduction, we note that the above ^"^x"*" ^'^'^^i^il^^s determining 

m~± and Arn^ (and, of course, detecting the charginos in the first place) imply a focus on events with 
the most energetic photons. Unless m~± is quite near y/s/2, the typical photon-tag energy and transverse 

Xi 

momentum will be large enough that two-photon (and certainly other) backgrounds to our 77r+7r~ 1^ final 
state will surely be negligible even if Arriy. is large enough that most events do not have an observable HIP 
for one of the vr's. In contrast, the techniques discussed in the following rely on using the full range of 
photon tag energies and transverse momenta, and could be compromised if (contrary to our expectation) 
there is significant background at lower E^,p7p values when Arrij^ is such that HIP(s) are not observable. 
We now compare the above Myy ~ 2m~± results to those obtained by employing the full distribution 

Xi 

in [M^^jE"*], including not only the location of the end-points in the E* spectra as a function of M^^, 
but also the full shape of the distribution. In particular, one could hope to make use of the large number 
of events at large M^^. We now describe two procedures for estimating the resulting errors on m~± and 
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AttIj^. In the first, we consider only the one-dimensional distribution in E* as a function of ^n-ji and Arriy.. 
We performed our study for rn^± values centered on 175 GeV and Am^ values centered on 200 MeV and 
employed bins of size 10 MeV (i.e. chosen to be somewhat larger than the resolution for E*). For two 
choices {A and B) of [?7^~±, Am^], we found the maximum difference -Dmax! based on the 1-D Kolmogorov- 
Smirnov (KS) |2^ procedure, between the two cumulative distribution functions as a function of E* and 
computed the probability, = P{D > Dmax), that the observed (or more precisely for theoretically 

computed event rates, the expected) value of Dmax is inconsistent with the two distributions having come 
from the same parent distribution rather than from two different distributions. The values of Q^g as a 
function of ^'m~± = m~±{A) — m~±{B) and 5Am^ = Am^(A) — Am-j^(S) are given in Fig. |^ for choices of 
m~± and Arriy in the ranges 173 < m~± < 177 GeV and 198 < Amy < 203 MeV. Due to the correlation 

Xi Xi 

apparent from Eq. (|5|), which shows that a change in Am^ can be compensated by a corresponding change 
in m~± without greatly affecting the E* distribution, we find that the more direct determination of m~± 

Al Al 

from the threshold region is crucial, especially for L = 50 fb^^. For L = 50 fb~^ (L = 1 ab~^), a 
determination of m~± to within 1 GeV (0.2 GeV) will allow a determination of Arrij^ to within ±5 MeV 
(±1 MeV) at the Q]^s = 0.3 (roughly la exclusion) level. 

A generalized KS test can also be applied to the full two-dimensional distribution in the [M^^,E*] 
plane. (As above, we employ bin sizes of 10 MeV for E*. For My-y., we employ bins of size 5 GeV, again 
somewhat larger than expected resolution error.) Results for the resulting Q'jpg values as a function of 
6m~± and SAra^ are given in the bottom two windows of Fig. ^. For L = 50 fb~^, we see that the statistics 
are somewhat marginal for use of the 2D technique, but a ±5 MeV determination of Am^ is possible if 
m~± is known to within 1 GeV. For L = 1 ab~^ the 2D KS test is superior to both the ID KS test based 

Xi 

on binning only in E* and the near-threshold technique. Using only the 2D KS test, we find that both 
m~± and Am^ can be measured with good precision. In particular, the Q'jps > 0.3 region is confined 
to 6m~± < 1 GeV, —1 < 6Amy < 2 MeV. If we combine the 2D KS analysis with the threshold 
determination of w-^i , we do even better. In particular, the threshold determination of n^^± to within 
0.2 GeV in the L = 1 ah~^ case, implies a reduced range of — 1 < 6 Amy < 1 MeV at the Q'jpg = 0.3 level. 

Let us now return to the impact parameter distribution for the vr^ in the final state and the extent 
to which it can be used to determine Am^, cr and underlying SUSY parameters. For Am^ ^ 800 MeV, 
the 7-|-IIIP(s) signature typically has substantial rate. We focus on the distribution of vr^ transverse 
impact parameters, 6t, ^ using the central reference case of it^^± = 175 GeV, Am^ = 200 MeV and 
cr = 23.81 mm (the nominal value from Table ||, based on Ref. for the given Am^ value). We generate 
e~^e~ — > jxtxi ~^ 77r+7r~^ events at ^/s = 600 GeV, using the minimal pj. > 10 GeV cut, for given input 

values of m~±, Amy and cr. We decay the xt according to the chosen cr, accounting for time-dilation 

Xi ^ 

for a given Xi^ velocity on an event-by-event basis. We compute the 6^ of the final state tt's if they pass 

through the toy-model vertex detector as specified in |11], which includes a first (LOO-like) layer at 1.6 cm. 

If bx is > 5ab, where ab is the pj^-dependent impact parameter resolution given in ||Tl|, then we enter 6^ 

'^Even though the exact location of the event vertex using the 7 trigger photon may be somewhat uncertain, bx can be 
rehably measured as the distance in the transverse plane by which the projected p,r unit vector misses the z axis. 
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Figure 6: We show different levels of Q]^g and Q^g computed from the E* and [M^^,E*] distributions 
(respectively) for two different choices, A and B, of ^n-^i and Am^, with 173 < rn~± < 177 GeV and 
198 < Amy < 203 MeV, as a function of 5m~± = m~±{A) - m~±(B) in GeV (taken to be > by 



convention) and 6Am^ = Amy.{A) — Am^{B) in MeV. Results shown are for L = 50 fb^^ and L = 1 ab^^. 
The symbols x, o, and □ correspond to Qks < 0.1, 0.1 < Qks < 0.3, and 0.3 < Qks < 0-68, respectively. 
Overlap of symbols means that different results are obtained for different choices of m~±{A) within the 

Xi 

above range. Tiny dots are points that cannot be distinguished, while the blank locations were not sampled 
in our analysis. 



18 



2250 



-^2000 
I 1750 

ID 

^ 1500 

0) 
Q. 

^jn 1250 



0) 1000 

Q. 
V) 

o 750 
> 

^ 500 

ID 

"° 250 







1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

- +1 


- 
■ 


: + + 


_ 
■ 


t 


- 


r + 




-_ + 


- 


+ 




T 


_ 


+ 


- 


+ 




-)- 






- 


; X 




















1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



01 23456789 10 

by (mm) 



Figure 7: For y/s = 600 GeV, m~± = 175 GeV, Amc = 200 MeV and cr = 23.81 mm, we plot the 6t 
distribution for the final vr^ in e~^e~ — > 7X^x1 after the cuts described in the text. The errors shown are 
those for L = 1 ab~^ of integrated luminosity. 



in the appropriate impact parameter distribution bin. (There can be up to two entries per event.) The 
resulting distribution for ?7ij± = 175 GeV, Arriy- = 200 MeV and cr = 23.81 mm is plotted in Fig. 0. The 
bulk of the bx distribution lies in the 6t < 10 mm range. This is also the portion of the range for which 
fluctuations will be under control for bins of reasonable size. We have chosen a bin size of 0.25 mm and, 
in what follows, we will consider only the bins with < 10 mm. Because of the complicated kinematics 
(including the trigger photon), the limited size of the vertex detector, and the 6^ ^ ^<^b{PT) requirement, 
the shape of the bx distribution is not a simple exponential, and, in particular, is cutoff at low 6^. Assuming 
all these effects can be adequately modeled and studied via Monte Carlo, the 6t spectrum can provide 
considerable information regarding the underlying parameters. Indeed, the spectrum has a complicated 
dependence on 'ni~± (as it affects the amount of vr boost), Am,^ (as it affects the amount of momentum 
available for motion transverse to the Xi direction) and cr. Before the bx > 5o";,(pt) requirement, the 
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bx distribution shifts to larger bx with increasing cr, decreasing Am^ and decreasing 'm~±. However, the 
br > 5(76 (pt) requirement removes more entries at low br for small Am^ than for large Am^, with the 
result that larger Am^j^ actually has more weight at small bx than smaller Am^. We now summarize the 
ability to use the br distribution to determine cr and Am^ for our detector model. We reemphasize that 
use of this technique requires that the acceptance and resolution of the vertex detector and the influence 
of the precise cuts made all be well understood. 

Vs = 600 GeV, Ax^ from by Distribution Shape 
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Figure 8: For 'm~± = 175 GcV, y/s = 600 GeV and central values of Am^ = 200 McV and cr = 23.81 mm, 
we plot the Ax^ computed as a function of cr — 23.81 mm for Am^ = 199, 200 and 201 MeV. We have 
assumed L = 1 ab~^, and have chosen the relative normalization of the alternative distribution for each 
alternative parameter set so as to minimize A^^. 

To be quantitative, we have proceeded as follows. We have chosen the central values of Am^ = 200 MeV, 
cr = 23.81 mm and m~± = 175 GeV. We then consider various neighboring values and compute the Ax^ 
between the bT distribution for any given set of [Am^ , cr, '"t.~± ] values and the central values using bins 
of size 0.25 mm in the region bx < 10 mm and choosing the relative normalization that minimizes Ax^ 
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(i.e. we rely on shape differences only). The distribution is insensitive to changes in m~± within the 
error of ±0.2 GeV from the kinematic distribution techniques. Thus, it is most relevant to assess our 
ability to determine cr and Amy from the brp distribution. Holding m~± = 175 GeV fixed, the Ay^ 

^ Xi 

values as a function of cr — 23.81 mm are plotted for Am^ = 199, 200 and 201 MeV in Fig. ^. For 
two parameters, the la (68.27% CL) and 2a (95.45% CL) Ax^ values are 2.30 and 4.61. We see that, 
assuming L = 1 ab~^, a Ax^ value below 4.61 can only be achieved for Am^ values within < 1 MeV from 
Amy. = 200 MeV. Further, it is apparent that cr will be determined to better than 5ct ~ 0.4 mm at the 
la level, corresponding roughly to a ±2% measurement of cr. 

It is perhaps useful to see what features of the bj- distribution allow such precise parameter determi- 
nation. To this end, we plot in Fig. |^ the differences, AA^, in the number of entries in a given bT bin 
between the distributions for the central [Am-j^ = 200 MeV, cr = 23.81 mm] parameter choice and those 
for [Am^^ = 201 MeV,cr = 24.11 mm] and [Auiy- = 199 MeV,cr = 23.41 mm]. The cr choices for these 
latter two Aniy. values are those which minimize the Ax^, as shown in Fig. ^. The AA^ values are those 
obtained after choosing the relative normalizations of the alternative parameter distributions so as to min- 
imize the Ax^. Also shown are the errors for a total luminosity of L = 1 ab~^. We observe a systematic 
variation, especially at low b^jp, of the Am^ = 201 MeV and Am^ = 199 MeV distributions relative to that 
for Am^ = 200 MeV. There are a large number of bins, each with Ax^ ^ 0.5, which combine to give the 
significant overall Ax^ plotted in Fig. ^. (We have checked that these differences between the distributions 
are not an artifact of the Monte Carlo simulation statistics.) Note for instance the depletion in the smallest 
6t bin for Aniy. = 201 MeV, mentioned earlier, followed rapidly by an excess in the next set of bx bins. 
We hope that these distributions make it apparent that both can be distinguished at a good level from 
that for [Am^ = 200 MeV, cr = 23.81 mm] and that the distributions for the two alternative parameter 
choices can also be clearly distinguished from one another. 

In the above analysis, we have implicitly assumed that there are no contaminating background events. 
After our minimal pj^ > 10 GeV tag cut, backgrounds (including the two-photon backgrounds) are unlikely 
to yield events having significant impact parameters.^ 

A measurement of cr can be converted to a joint constraint on Am^ and the underlying Mi , M2 , /i, tan /? 
parameters. Consider Am^ < 700 MeV for which the only modes of any importance are xt ~^ ^^z^Xi 
Xi — > TT^Xi- For the latter, we have 



m~± — m~o 

Xi 



Xi 

[oh - [(m|± - m|o)' - ml (m~± + m~o 



(7) 



where the O^l^ describe the W^xfxi coupling. 



On = -^NiiVi2 + ATiaVii , Of^ = +1=NM2 + NMi • (8) 



^If cr is of the order of 100 /im, then a precise accounting for the e'^e — > e'^e 7t^t — > yn^Tv Ff background might be 
necessary. 
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Figure 9: For m~± = 175 GeV and ^/s = 600 GeV, we plot the differences, AA^, (of the number of entries 
in a given bin) between the bx spectra for the central values of [Aniy. = 200 MeV, cr = 23.81 mm] and 
those for [Am^ = 201 MeV, cr = 24.11 mm] and [Am,^ = 199 MeV,cr = 23.41 mm]. These differences 
plotted are those found after choosing the normalizations of the distributions for the latter two alternative 
parameter choices so as to minimize Ax^. Also shown are the errors for a total luminosity of L = 1 ab~^. 
We have artificially shifted the bx values at which the points are plotted so that the two different cases 
can be distinguished. 



Current experiment gives = 92.42ib0.3 MeV, which yields a ~ 0.7% error for this width. More important 
is the parameter dependence of . Naively, these are ~ 1 for the wino-like LSP scenario. However, 
their dependence on the underlying SUSY parameters Mi, M2, /U, tan /3 is significant. For example, for 
m~± ~ 175 GeV, they vary by about 5% as |^| increases in the range 500 GeV to 1 TeV, with much 
more dramatic variation as approaches 300 GeV. At large |^| > 500 GeV the 0-^{ also vary by about 
(1 — 2)% as tan/3 ranges from small to large values. At smaller the tan/3 variation is much more 
dramatic. The widths for the Xi ~^ ^^'^Xi (^ = /x) modes are also proportional to combinations of the 
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Oii^. If the Oii^ are held fixed at the Ref. central values, Eq. (^) implies that a shift in Am-j^ by 
1 MeV corresponds (in the vicinity of Aniy. 200 MeV) to a shift in cr by about 2.5%. However, there 

L Ft 

are many parameter choices that yield the same Am^ with Oi{ values that differ by more than this. 
Given these uncertainties and the ~ 2% experimental error for the cr determination, it will be difficult to 
constrain Arriy. more accurately than via the experimentally direct kinematic distribution and shape 
fit techniques. Ultimately, the reverse strategy might prove useful. That is, determine Am^^ directly from 
the kinematic distributions and bx distribution shape and cr from the bx distribution shape and then use 
these values to constrain the Mi, M2, /i, tan/3 parameters. Note that it would be crucial to include the 
one-loop corrections to Am^ in this process. 

5. Determining M2, u and tan/3. Let us suppose that we have accurate determinations of m~± 

Xi 

and Arriy. using the kinematic distribution techniques described above. We will then wish to extract the 
underlying SUSY parameters. We study our ability to do so for the specific case of m~± = 175 GeV and 

Xi 

Am^ = 200 MeV. This particular case is motivated in the context of {Sgs = 0-II)/AMSB one-loop 
boundary conditions: (a) m~± < 200 GeV is preferred if we are to avoid the extreme fine-tuning for larger 

Xi 

m~± implied by the large value (~ 10) predicted for M3/M2; and (b) Am^ < 200 — 500 MeV is preferred 
because of the large value predicted for M1/M2. We again remind the reader that, for such Am^, the 
IXiXi events we employ will be background free since at least one of the final tt's will have an observable 
HIP. 

The SUSY parameters entering into the chargino sector at tree level are M2, fJ- and tan /5. The measured 
value of m~± provides only one constraint on these three parameters. The projection of this constraint 

Xi 

on the standard [/x, M2] plane gives the two regions indicated in Fig. [l^ by the outer solid lines. In this 
paper, we restrict ourselves to > 300 GeV, since for lower the xt charginos would also be produced 
and the nature of the analysis and parameter extraction procedures would change dramatically. Typical 
radiative electroweak symmetry breaking favors large |/i| values. 

The goal is to use the 'yxtxi production cross section and kinematic dependencies to obtain additional 
constraints on the three parameters. Sensitivity to these parameters arises entirely through the coupling 
of the Z to XiXi- 

^j>^(0[iPL + 0[fP^) , with Olf = -c^ + , 0[f = -c'^ + lu!2, (9) 

where cw = cos^vy etc., and V12 and U12 are elements of the matrices that diagonalize the chargino mass 
matrix. (In our work, we neglect possible CP-violating phases in the chargino sector.) Their squares can 
be written in the form = ^(1 — cos 20/,) and U12 = ^(1 — cos 20/j) with 

M| - ^2 - 2m^ cos 2(3 



cos 20/, = - 

M|-/x2 + 2m2, cos 2/3 

cos 2(f>R = , . (10) 



(M| + /x2 + 2m^)2 - 4(M2/i - m^^ sin 2/3)2 

M| - + 2m^ cos 2/3 
Ml + /x2 + 2m2^)2 - 4(M2/i - m^, sin 2/3)2 
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Figure 10: For m~± = 175 GeV, the outer solid lines show the region in the [/U, M2] plane consistent 
with 1 < tan/3 < 100 and |/i| > 150 GeV. Also shown are the values of M2] sampled for tan/? = 
2, 5, 10, 20, 35, 50 in determining our fits to the e"^e~ — > jxtxi cross sections. 



Asymptotic expressions for Vu and U12, valid for IM2 ± /i| mz are ||2l| (using c/j = cos/3 etc.): 



M|^7 ' ■ ^ ^ 

Clearly V12 and ^'^^ small compared to when is large. Further, UI2 3> V12 is typical for tan/3 > 2 
when ^ M2. In the scenarios we consider, the Xi ^'^^ indeed highly wino-like and, thus, the problem 
is to pick out the small Vy2 and UI2 corrections to the dominant —c^/ term in the Z coupling. 
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To understand how to proceed, it is useful to first briefly review results for e~^e~ 
photon tag. Following [^2[, we write 

o2 



Xi'Xi without a 



XiXi 



e 
s 



a,l3=R,L 



(12) 



where 



Qll{s) 
Qlr{s) 
Qrl{s) 
Qrr{s) 



1 + 



1 + 



1 + 



1 + 



Dz{s) 
Dzjs) 
Dz{s) 
Dz{s) 



2 

{si. - \) 



12 



4s2 ^ 



~'av + 2^12 



^2 -U -r/2 
'-■w ^ 2 12 



In Eq. dl 

Q 



It is crucial to observe that Q rl 



^V?2 and 



RR 



2ct 



Dz{s) = s/is-ruz) and D~{s,t) = s/{t-m~ 

at large s, implying that a right-handed polarized e~ beam will provide a very direct probe 
of and [7^2, but at the sacrifice of a very suppressed cross section. Pure e]^ also has the advantage of 
eliminating the 'background' from the exchange diagram. Since this diagram will significantly suppress 
the unpolarized cross section for < 1 TeV (i.e. even for masses for which we will not be able to directly 
detect the t'e and measure its mass), the normalization of the unpolarized cross section cannot be calculated 
reliably enough to use as an ingredient in extracting the parameters of the model. The importance of the 
polarized measurements, for which this is not a problem, implies that maximal integrated luminosity will 
thus be of paramount importance, whether one is looking at an mSUGRA-type SUSY breaking scenario 
or, as we shall discuss in detail shortly, the degenerate m~± ~ m~o scenario. In terms of the Qafl's, the 
e+e~ — > xfxi amplitude squared is 



OC (-2p~--Pe-p~+-Pe+s(QLi? + Q|L)-2p~7-Pe+Py^ 

(13) 



where p~- ■ p, 

Xi 

Eh 



PxfP- 



PCOSO), p~- -Pe 



PxfP- 



1 



m-^/Er and 

■y/s/2. By measuring the magnitude of the cross section, especially for a right-handed polarized 
beam, and its dependence on cos 9, one can hope to determine the charges Q^p and, thence, the crucial 
and U12 matrix entries that probe /i, M2 and tan /?. The Xi are not directly observed (unless they have 
a substantial path length, as is possible in this model but not assumed in our analysis), but, by using the 
accurately measured values of m~± and Amc and the measured three- momenta of the vr^, the chargino 

Xi 

momenta can be reconstructed up to a two-fold ambiguity. The analogue of this reconstruction for the 
IXiXi final state is discussed in the Appendix. We note that \A\'^ and, thence, the integrated cross section 
for any particular machine energy and set of cuts depend bi-quadratically on V12 and [7^2 ■ 

These same 'charges', Qap, appear in the expression for the e+e" ixtxi amplitude, which contains 
two initial state radiation diagrams summed over QapiQ^ = -^xx^'® final state radiation diagrams 

summed over Qapio^ = s)'s, where q is the four-momentum carried by the virtual Z or 7. The resulting 
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Vs = 600 GeV, ni^+ = 175 GeV 





45 








44 


unpol 


43 


b 






42 




41 




40 




39 




38 



^ 0.32 

O 

a. 
c 

< 0,315 
0.31 
0.305 
0.3 
0.295 
0.29 



' ^^^^ 



X 
X 
X 



- X ianlS=2 
■ ianl3=50 



•1000 -500 



^^^^^^ i 



X I 



b 



500 1000 

M (GeV) 







; X 




— X 


x 


X 

— XB 
X* 


*^ 

B 


- X >^ 


bX 

X 


:■ X j B 

- ■•?<^ 

-1 "b I 


>*B^ 

B 'y^^m 

■ B X 



< 



•1000 -500 



500 1000 

/J, (GeV) 



0.2 
0.175 
0.15 
0.125 
0.1 
0.075 
0.05 
0.025 


0.8 
0.6 
0.4 
0.2 


-0.2 



X 

B 

X 



•1000 -500 



500 1000 

^ (GeV) 



- \ 

X X 


x\ 

X 

B 


■ 

X ■ 


B 
B 

X 


_B _ 

X B 


■ 


X B 

■ B 


X 

bX 

X 

|X , 



1000 -500 



500 1000 

/X (GeV) 



Figure 11: For m~± = 175 GeV, we plot a and A as a function of ri for tan/3 = 2 and 50. Results are given 
for an unpolarized e~ beam and a purely beam at ^/s = 600 GeV. A mass splitting of Am^ = 200 MeV 
was employed in generating the vr^ momenta. 

form for \A\'^ has a complicated dependence on many kinematical variables. (In the Appendix, we give the 
structure of \A\'^ for the dominant initial state radiation diagrams. An expression in a different formalism 
for the full is given in |]l5[.) In the soft-photon limit, |^(e"'"e~ —>■ 7Xi'xr)P directly proportional 
to \A{e~^e~ —>■ Xi'xr)!^' but this limit does not allow an understanding of all of the features of the 7X1" XF 
cross section. However, it remains true that right handed electron polarization produces a very suppressed 
cross section with high sensitivity to and and no sensitivity to the probably unknown z/e mass. 
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We now outline our precise procedures. For each event, the observed 7 defines a chargino-pair, i.e. 
M^^, center of mass system. In this c.m.s., we have fohowed the Colhns-Soper procedure [23| of defining 
a z-axis by boosting along the 7 direction to the My.^ c.m.s., determining the unit vectors u^± in the 

three momentum directions after the boost and defining z = {u^+ — u^-)/{2 — 2u^+ • u^-)- For the 
vast majority of events, z is very closely aligned with the z axis in the laboratory frame. We then 
determine the xt momenta in this frame (using the observed vr^ momenta and the known values of m~± 
and m~o = m~± — Amy) following the procedure of the Appendix and compute the angle, 9, of p~~- with 
respect to z. {p~+ = —p~- in the M^^ c.m.s.) Since there is a two-fold ambiguity in the reconstruction, 
we bin events twice corresponding to the two possible reconstructed xt three-momenta. 

We have restricted our analysis to the absolute magnitude of the cross section, a, and the asymmetry 

ufcos ^ > 0) - cjfcos ^ < 0) ^, 
A = — (14) 

cr(cos >0)+ a{cos 9 < 0) 

where the cross sections are obtained by integrating over all of phase space consistent with the specified 

sign of of cos^ and our basic pj, > 10 GeV and 10° < 0^ < 170° photon-tag cuts. We illustrate the 

dependence of a and A on fi for tan (3 = 2 and 50 in Fig. |ll|. (M2 is determined as a function of /i at a 

given tan/? by the fixed value of m~± = 175 GeV as shown in Fig. plH.) Results are given for unpolarized 

Xi ' — ' 

beams and for the case of a purely right-handed polarized electron beam. For L = 1 ab~^, typical la errors 
for a and A are of order 0.5% and 1% (15% and 50%) for the unpolarized (polarized) cases, respectively. 
Relative to the above-quoted typical errors, there is significant variation with respect to /u of both a and 
A, but different tan /? values will be difficult to separate since a large change in tan /? can be compensated 
by a small shift in fi that will leave the four observables more or less unchanged. The variations with fj, are 
especially dramatic for iTpoiarizcd and Apoiarizedj but the large statistical errors even for L = 1 ab~^ render 
these variations no more (or less) useful than those for the corresponding observables in the unpolarized 
beam case. Since the polarized and unpolarized cross sections and asymmetries contain different weightings 
of U12 and V12, accumulating significant luminosity in both modes is useful for maximizing our ability to 
determine and Vi2- 
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-7.08801 
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at 


0.0111671 


-0.323629 


-0.124093 


2.91672 


0.637143 


1.26579 


ap 


0.00771917 


-0.103494 


-0.169943 


0.59423 


1.27397 


0.784927 



Table 6: Tabulation of fitted parameters a, b, c, d, e, f (in units of fb's) of Eq. ( p^ ) for each of the four cross 
sections employed in our analysis. Results are for ^/s = 600 GeV and are those obtained after Monte Carlo 
integration incorporating cuts and the chargino reconstruction algorithm described. 



To simplify our analysis, we scanned a large selection of // values in the ranges [—1000 GeV, —300 GeV] 
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and [300 GeV, 1000 GeV] at a series of tan/3 values (tan/3 = 2,5,10,20,35,50)0 and fit the results for 

1) = 0-unpolarized(cos6' > 0), 2) CJ^ = CTunpolarizcd (cOS < 0), 3) (J^ = CTpolarized (cOS 6* > 0) and 4) 

dp = cTpoiarized (cos < 0) to the bi-quadratic form (required by the amplitude structure) 

ai = a + hUl^ + cVil + d{Ul^f + e{V^^f + fUf^V^\ . (15) 

Results for a — / for ^/s = 600 GeV, m~± =175 GeV and Am^ = 200 MeV (the latter possibly affecting 
reconstruction of the angle of the chargino) are given in Table ^. These fits give the correct cross section 
with an accuracy of better than 0.7% for all the sampled points. For ap and cxp, note the very small 
constant terms and the much larger quartic term coefficients {d,e,f). As already discussed, this is to be 
expected and follows from the wino-like nature of the xf . 

We now employ these cross section fits to study the accuracy with which we can expect to determine 
M2, n and tan (3. To illustrate, we give in Fig. |l|, for ^/s = 600 GeV, contours of Ax^ = 2.30 (68% CL for 
two independent parameters) and Ax^ = 13 (99% CL), using reference models of [tan/3,/i] = [5,375 GeV] 
and [5, 600 GeV]. For this figure, we employ the following set of observables: the forward/backward asym- 
metry for unpolarized beams; the absolute cross section for polarized beams; and the forward/backward 
asymmetry for polarized beams.0 As explained later, the absolute unpolarized cross section is likely to have 
systematic theoretical uncertainties that are much larger than the statistical errors, and is not employed. 
We observe that, at the 68% CL, M2 and \fj,\ can be determined to ~ 8% and ~ 16%, respectively, for 
1^1 ~ 400 GeV and to ~ 8% and ~ 40%, respectively, for |;u| ~ 600 GeV; tan/3 is essentially undetermined. 
The 68% CL limits on /i and M2 for a wider selection of input /i values and input tan f3 values of 5 and 20 
are shown in Fig. Results for other tan (3 values are quite similar. Unless tan /3 has been determined by 
other data, we must combine all such graphs and take the outer errors. Clearly, at the highest |^| values, 
errors for the fj, and M2 determinations become quite large, and significant uncertainty in M2 develops as 
fi falls below 300 GeV. 

The unpolarized cross section can provide useful constraints on SUSY parameters only if systematic 
errors are substantially smaller than the typical 5%-10% variation of Uunpoiarized as a function of parameters; 



see Fig. |11|. Systematic experimental errors at an electron collider should be quite small, probably even 
smaller than the ~ 0.5% statistical errors. However, theoretical uncertainties in (Tunpoiarized must also be 
well below the 5%-10% level. Even if we assume that the t'e contribution can be absolutely normalized (for 
example, using z^e observations from the LHC) higher-order electroweak corrections to the e"*'e~ 'yxtxi 
cross section will need to be computed. In addition, one will be implicitly assuming that the supersymmetric 
ZxiXi coupling strength is indeed precisely that predicted by assuming strict supersymmetry for tree-level 
couplings. This will be difficult to independently check. Finally, to employ the absolute normalization of the 
unpolarized cross section one would need to rely on the modeling employed for computing the Xi ~^ ^^Xi 
branching ratio at the ^ 5% level. The computation for the xf ~^ ^^^Xi xt ~^ ''^^Xi widths was 



^The value of M2 at each point is determined by m~± . 

^1 



''We compute the error in A= " ' ~" as SA = J where N ^ N+ + N 
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Figure 12: For m~± = 175 GeV, we plot the contours for Ax^ = 2.3 (68% CL for two independent 
variables) and Ax^ = 13 (99% CL) in the [/U,M2] and [/i,tan/?] parameter spaces using the two reference 
cases of [tan /?, ^] = [5, 375 GeV] and [5, 600 GeV]. We have assumed running at -y/s = 600 GeV sufficient to 
accumulate integrated luminosities of L = 1 ab~^ for both unpolarized beams and for a purely right-handed 
electron beam. The Ax^ has been computed using data for (Tpolarizedi ^unpolarized and ^polarized? '^unpolarized 
was not included. A mass splitting of Aniy^ = 200 MeV was used in generating the pion momenta. 

reviewed in the previous section. While some of the theoretical uncertainties in the widths cancel when 
computing the branching ratios, the latter will still have dependence on the Mi, M2, /i, tan /? parameters 
at the (2 — 5)% level for Am^ ^ 700 MeV (much larger as the multi-pion decay modes of the Xi enter at 
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Figure 13: For m~± = 175 GeV, we plot the upper and lower Ax^ = 2.3 (68% CL) limits for and 
M2 as a function of the input value of fi, fiQ, for input values of tan/3 = 5 and tan/3 = 20. The dashed 
lines indicate either the input values of fiQ or the corresponding M2(^o) values for the input tan/3, given 
m~± = 175 GeV. For each /io, the 68% CL boundaries for choices of fi with the same (opposite) sign of fi 
are shown by solid (dotted) lines. Points not plotted at small /i ~ 300 GeV are inconsistent with the given 



input m~± 



175 GeV. Data inputs are as specified in Fig. 12 



somewhat larger Am^). For this paper, we have assumed that the absolute normalization of the unpolarized 
absolute cross section is not computable to the 5% level, and have determined the Ax^ values using only 
^unpolarized, ^-polarized and ^polarized- If we Were able to reliably use dunpoiarized , the resulting plots analogous 



30 



to Fig. ^ would show a moderate increase in our ability to determine M2 and /x, but tan {3 would remain 
essentially undetermined. 

We have also examined whether it is useful to include data from several ^/s choices for determining 
/X, M2,tan/3. We found that splitting the total luminosity between ^/s = 450 GeV and 600 GeV actually 
causes a mild decrease in our ability to determine the parameters (assuming no large variation of instan- 
taneous luminosity with y/s) . This remains true even if we run at both energies and use as an observable 
the ratio o"unpoiarizod(\/^ = 450 GeV)/(Tunpoiarized(\/s = 600 GeV) (which would be less affected by system- 
atic theoretical uncertainties as compared to the individual absolute o"unpoiarizcd normalizations). If we are 
conservative and rely only on polarized data (L = 2 ab~^ at ^/s = 600 GeV), the M2 and /i accuracies are 



almost the same as shown in Fig. 12 



Additional information can, in principal, be gleaned from the Xi^ ~^ Xi"^^ decays. However, this is 
quite difficult in practice. The issue is whether one can use correlations involving the pions in the final 
state. Correlations between the pions themselves or between the pions and the reconstructed momenta 
for the Xi^ will be present to the extent that information regarding the direction of the spins of the 
decaying Xi^ influences the vr^. Consider Xi decay. For a given s~ choice, one finds 

l-M(xr - vr-x?)P oc 1 + ^-^x-Pj--^- (16) 

'^Xi 

after dropping terms higher order in the small ratios /^m^/m~± and PTt-/m~±. Since is of order 

Am^, the dependence on p^- ■ s~ is obscured by the much larger leading term. For this reason, we have 
not attempted to employ these correlations. 

6. Discussion and Conclusions. We have considered the detection at a linear e~^e~ collider of the 
lightest chargino in a model in which it is closely degenerate with a wino-like neutralino. The phenomenol- 
ogy depends critically upon the mass splitting Am^^ = '"^~± — ^xl' '^^^ ^"^x ^ ""^t' ^^^^ observe 
long-lived heavily ionizing Xi^ tracks. For Am^ ^ 2 GeV, the usual mSUGRA jets -|- missing energy 
signal will be dominant. Both can be detected for m~± essentially all the way up to the threshold at 

Xi 

\/^/2 assuming modest integrated luminosity of L = 50 fb^^. For m~± not near threshold, much less 
luminosity would be adequate. For the range 200 MeV < Am^ < 1 GeV, favored in typical models, the 
7 -|- vr+vr^ -|- Ip final state will be the crucial discovery mode. Then, so long as soft > 200 MeV pions are 
visible in the detector ,Q discovery in the ■y + 7T~^TT~ + Ip mode will be possible up to 'ni~± ~ {y/s —pj,)/2 for 
L = 50 fb^^. In particular, we argued that simple cuts requiring pj. > 10 GeV, only two soft/central pions, 
large (mainly invisible) mass recoiling opposite the trigger 7, and no small angle {9 > 1°) e~ or (as 
would be present for the most obvious two-photon backgrounds) have an excellent chance of suppressing 
backgrounds to a negligible level. An overview of the discovery modes relevant as a function of location in 
the full ['m~± , Am^] parameter space appeared in Fig. 



We then turned to a more detailed study of scenarios that would arise in the 5gs = O-II string model 



^If the soft pions are not visible, 7X1^X1 production yields a 7+ f] final state with substantial background from 'y + Z{^ uu), 

and discovery reach will be more limited and the error for m~± , summarized shortly, much larger; Am^ will not be measurable. 

^1 
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and in the AMSB model. These models predict the same low-energy ratios for the gaugino masses, M2 : 
Ml : Ms ~ 0.3 : 1 : 3 (tree-level), and generally require substantial |/x| for automatic electroweak symmetry 
breaking. In these models, the smallness of M2 compared to the other mass parameters implies that Am^ 
will almost certainly lie in the m,r < Am^ < 800 MeV range (after including radiative corrections) for 
which the only appropriate mode for discovery and study is e^e~ — 'yXiXi ~^ l'^^'^~XiXi- The large 
ratio of M3/M2 implies that fine tuning for these models will be quite extreme unless the chargino mass 
lies in the < 200 GeV range. For m~± and Am^ in these ranges we studied the accuracy with which 
m~± and Am^ can be measured using kinematic distributions of the pions in the final state. For the 
particularly typical choices of m~± = 175 GeV and Am^ = 200 MeV, we studied the accuracy with which 
the Am^ and cr could be determined from the transverse impact parameter distribution for the soft pion 
tracks. Finally, for m~± = 175 GeV and Amy = 200 MeV, we determined the precision with which the 
fundamental underlying SUSY parameters M2, n and tan/3 could be determined from cross sections and 
asymmetries. We summarize the results of these studies below. 

The errors estimated for m~± and Am^ assume that the electron sneutrino is not so light as to sub- 
stantially suppress the unpolarized JXiXi cross section. Error estimates for m~± and Am^ obtained by 
studying kinematic distributions of the final pion decay products were given for two procedures, assuming 
L = 1 ah~^ of accumulated data. The first employs events in which the XiXi invariant mass is very close 
to 2m~± (i.e. events with very large pj,). After accumulating L = I ab^^ with unpolarized beams, we 
find that the location of the threshold for the recoil mass opposite the 7 trigger (i.e. the invariant mass 
of the XiXi P&ii' system), denoted M-^-^, gives a very accurate measurement of ^^i-i (typically < 0.5% 
error), and that the average soft vr energy, E*, in the center of mass for My.y. near threshold gives a 
very accurate measurement of Arn^ (typically 1%). These procedures for detecting 'yXiXi ^iid measuring 
m~± and Aniy. using only events with the highest photon energies, such that M^^. ~ 2m~±, guarantee that 
these goals can be achieved even if there is significant background (contrary to our expectation) at lower 
E-f and pj^ for the photon tag. Only if "i-^± is close to would we have to rely on the background 

being small for pj, as small as 10 GeV. However, if we are correct that this background remains small for 
such a low for the trigger, slightly better accuracies for rn~± and Am^ can be achieved (for all, hnt 
especially moderate, '>n~± values) by employing the full [E*,M^-^] distribution of the events and applying 
statistical tests for discriminating between the distributions obtained for different [?^^±5 Anij^] choices. We 
emphasize that the low pj, events will certainly be background free if Am^^ < 600 MeV (as typical) since 
in almost all events one or both of the final tt's will have an observable high impact parameter (HIP). 

For the typical case of m~± = 175 GeV and Am^ = 200 MeV (for which essentially all events have a HIP 
TT and, hence, backgrounds are certain to be small even when including all events with pj, > 10 GeV), we 
found that for L = 1 ab~^ the value of Am^ can be determined to better than 1 MeV from the transverse 
impact parameter distribution of the final pions. Simultaneously, cr for the xf ^^^Xi decay can be 
determined to roughly 2%. (The impact distribution is less sensitive to m~±, which is best determined 
from the kinematic distributions.) 
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Finally, the underlying parameters M2 and ^ of supersymmetry can be most reliably determined from 
the magnitude of the "fXiXi right-handed polarized e~ beam cross section and the polarized beam and 
unpolarized beam angular asymmetries. As part of the procedure, one inputs the precise determination 
of m~± obtained by employing the threshold or statistical techniques. For the typical m~± = 175 GeV, 
Am^^ = 200 MeV case, we determined the statistical errors for these quantities assuming L = 1 ab~^ 
each in unpolarized and polarized running. We then computed the Ax^ for discriminating between a 
given input model choice of M2, fi,taji (3 (constrained to give the value of m~± = 175 GeV, presumed to 

Xi 

already be well-measured) and other possible choices consistent with the same ?ra~± value.0 We found 
the best parameter accuracies by accumulating luminosity only at the highest energy (y/s = 600 GeV 
in our study), distributed roughly equally between unpolarized beam running and pure running. We 
found la (68% CL) accuracies for M2 and of ±8% and ±16%, respectively, for |^| ~ 400 GeV, and 
±8% and ±40%, respectively, for |^| ~ 600 GeV. Errors for both M2 and \fi\ become uselessly large by 
1/^1 ~ 900 GeV. The sign of fi is not determined at the la level. However, the error for M2 decreases 
to about 5% for 350 GeV < < 750 GeV if the sign of fi is known from other input. Finally, tan/3 
is essentially undetermined at the la level. We have noted that for the small values of Am^ natural in 
the wino LSP scenarios, correlations involving the final soft vr's are negligible and do not aid in parameter 
determination. 

In the above, we did not employ the absolute normalization of the unpolarized cross section since it 
is sensitive to many theoretical uncertainties at the (5 — 10)% level, a level of uncertainty that is larger 
than the amount of variation with respect to the parameters of interest. One source of uncertainty is 
the unknown mass. By measuring the unpolarized 'yxtxi cross section at two energies, could be 
extracted with an accuracy determined by the other theoretical uncertainties in i^unpoiarized- However, if 
we expend luminosity for this purpose, the errors for the M2 and fj, determinations would increase. 

Of course, additional SUSY signals will emerge if some of the squarks, sleptons and/or sneutrinos are 
light enough (but still heavier than the Xi) that their production rates are substantial. In particular, 
leptonic signals from the decays [e.g. if^ ^^Xi or I'e ^^Xi ] would be present. Also, depending upon 
^/s and the mass splitting between the wino- like Xi and the bino-like X2 (which is large in many models), the 
suppressed X1X2 production channel might be detectable. Or, the radiative electroweak symmetry breaking 
scenario chosen by nature could be sufficiently unconventional that |^| is small enough for production of 
charged and neutral higgsino-like states to be detectable. We have chosen to emphasize the case (which 
is most likely in typical models) that none of these additional signals are present until ^/s substantially 
above 600 GeV is available. 

It is, of course, likely that the LHC will have been operating for a number of years prior to the 

contraction of the next e^e'^ collider. Because of the large available center of mass energy, it is probable 

that some of the heavier states mentioned above will be produced. However, it is not completely clear that 

they will be observed. The role of the LHC for the type of model being considered will be discussed in a 

^Values of < 300 GeV were not considered; for such values, xt production becomes possible and the strategy for 
determining SUSY parameters would change substantially. 
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later paper. 
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7. Appendix 

In this Appendix we present some details regarding the reconstruction of the momenta of all the particles 
for the process e~^e~ 'yxfxi ~^ 7^'''Xi'''"~Xi- (We use the notation to distinguish the Xi associated 
with the Xi decay; the Xi is, of course, its own antiparticle.) We also give the form of the jxtxi cross 
section resulting from the often dominant terms (a gauge invariant subset) in which the 7 is radiated only 
from an initial or e~ . 

The reconstruction of the final state momenta is performed in analogy with the techniques developed 
for e^e~ — > W~^W~ — > £^u£~T' [^]. We begin by first using the observed 7 four momentum and the 
known and e~ beam momenta to boost to the xtxi center of mass. The momenta in the expressions 
below are those defined in this frame. (We drop the * superscripts employed in the main text.) In the 
xtxi c.m.s., the 6 components of the observed 7r~^ and vr" three momenta, combined with the known m~o 
and m~± masses, can be used to solve for the three momenta of the Xi and Xi up to a two-fold ambiguity 
indicated by the ±d term below: 

p~o = ap^~ + bp^+ ± d{p^- xp^+ ), p~o = - {pzo + p^+ +Pn- ), (17) 

where a, b and d are given by:^ 

_ m|p;+p - np^- ■p^+ ^ -mp^- ■ p^+ +n\p^-\'^ 

\P-K- X Ptt+I K- X Ptt+I 

{Eb -E^^f - mK - \p^~ ? - b'^\p^+\'^ - 2abp^- ■ p^+ 
= , (18) 



where 



i (m|± - m|o - 2EbE^- + E^. + \p^-\^) , n = ^ (m|± - m|o - 2EhE^+ - 2p^- ■ p^+) . (19) 



m ^ 

2 

Using these formulae, we can, up to the two-fold ambiguity, reconstruct the xt ^^'^ Xi momenta. The 

construction only fails if p^~ x p^+ = or if < (the latter being possible as a result of momentum 

smearing in the detector). If < 0, setting d = usually gives a fairly accurate result for p~o and p^. It is 

perhaps useful to keep in mind approximations that follow from the fact that E^+ , E^- , \p^+ \ , \p^- \ ~ Am-j^ 

and m~+ — m~n ~ 2m~oAm^. For small Am^/m~n we have 
xf Xi Xi ^ ^' Xi 



E't-ml, 



m = EfjE^- — m~oAm^ , n = — EbE^+ + m-oAmy. , d = — — y. (20) 



\Pn- X P-n 



We correct some errors in the formulae of Ref. 124| 
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Further, p~+ = + p^- and p~~ = p~o + Pt,+ ~ p~o ■ Finally, we note that if > m~± Am^/m.^ , 

the TT directions are guaranteed to have positive dot product with the directions of their parent charginos. 

The utility of this reconstruction follows from an understanding of how the matrix element squared 
depends upon the final state momenta. As described in the text, correlations between the helicity of the Xi 
and the three-momentum of the decay pion are of order Am^/m~±. Thus, to a very good approximation 
the distribution of the decay pion in the xf rest frame is completely independent of the 0*, (p* rest frame 
decay angles for any fixed helicity of the xf. The final state only contains information encoded in the 
decay distributions of the xtiXi- Thus, aside from small terms of order Am^/m~±, the decay pions do 
not have any correlations other than those kinematically induced by the boosts in the directions of their 
parent charginos and the form of the cross section as a function of 9, the angle of the Xi iii the Xi'Xi 
center-of-mass. 

For the case of a photon radiated from the initial e"*" or e~ , the form of the invariant matrix element 
squared for e~^e~ — > ^XiXi is = (2pe+ 'P-yPe~ ' P-y)~^ I'Sp with 

|Bp OC -2p~- ■ Pe-P^+ ■ Pe-Pe+ ' P-yiQlL + QlR + QrL + Qrr) 
-2p~- • Pe-P^+ ■ Pe+S {QIr + Qrl) 

+2p~- • Pe~P^+ ■ Pe+ (Pe- ' P^QlR + Pe- ' P^QrL + Pe+ " P-yQlR + Pe+ ' P-yQ^l) 



+PZ- ■ Pe-SP~+ ■ P^iQLR + 



RL) 



• Pe-P^+ ■ PiPe- ■ PiiQhR + Qrl) 



2 , n,2 
LR ~ 

2 I 



-2p~- • Pe+P^+ ■ Pe-HQiL + Qrr 



+2p~- • Pe+P~+ ■ Pe- (Pe- ' P^/QlL + Pe- ' P-tQrR + Pe+ ' P^QIl + Pe+ ' PfQ^ 



RRl 



Xi 

-2p~- • Pe+P^+ ■ Pe+Pe- ' Pj{QlL + QlR + QrL + Qrr) 

+P^- ■ Pe+sp~+ ■ P-riQiL + Qrr) 

-2p~- ■ Pe+P^+ ■ P-yPe+ " P-yiQll + Qrr) 

+P^+ ■ Pe-sp~- ■ PyiQli + Qrr) 
-2p~+ • Pe-P^- ■ P-yPe- ■ P-riQiL + Qrr) 
+P^+ ■ Pe+sp~- ■ P-riQLR + Qrl) 

-2p~+ ■ Pe+P^- ■ PjPe+ ■ P-yiQLR + Qrl) 

-s'^m~±{QLLQLR + QrlQrr) 

+2s m~± (Pe- • P^QllQlR + Pe- ■ P-jQrlQrR + Pe+ ' P-/QllQlR + Pe+ ' P-^QrlQrr) 

-2m~± ((Pe- • P-t)'^QLLQLR + {Pe- ' P^i^QrlQrR + {Pe+ ' P^)^QllQlR + {Pe+ ' P-y)"^ Q RlQ Rr) ■ 

In the above, all Qap^ are evaluated at M"^^. Note that \A\'^ reduces to Eq. (13) in the — > limit after 
removing the photon radiation pole factors {p^+ • p-f){Pe- 'P'y)- This answer is gauge invariant on its own. 
Contributions to coming from diagrams in which the photon is radiated from one of the final Xi lines 
are often suppressed relative to the initial state radiation terms given above. All diagrams were included 
in our numerical calculations. In practice, there are delicate cancellations between large terms in Eq. (^T|), 
and also in the full expression for \A\'^ obtained after including all diagrams, that can lead to numerical 
inaccuracies. Thus, it is actually best to sum all the diagrams numerically at the matrix element level and 
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then square. This is the approach used in our numerical simulations. 
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